Magic numbers in exotic nuclei and spin-isospin properties of NN interaction 
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The magic numbers in exotic nuclei are discussed, and their novel origin is shown to be the spin- 
isospin dependent part of the nucleon-nucleon interaction in nuclei. The importance and robustness 
of this mechanism is shown in terms of meson exchange, G-matrix and QCD theories. In neutron- 
rich exotic nuclei, magic numbers such as N=8, 20, etc. can disappear, while N—6, 16, etc. arise, 
affecting the structure of lightest exotic nuclei to nucleosynthesis of heavy elements. 



PACS number(s): 21.30.Fe,21.60.Cs,27.20.+n,27.30.+t 



The magic number is the most fundamental quantity 
governing the nuclear structure. The nuclear shell model 
has been started by Mayer and Jensen by identifying the 
magic numbers and their origin |lj]. The study of nu- 
clear structure has been advanced on the basis of the 
shell structure associated with the magic numbers. This 
study, on the other hand, has been made predominantly 
for stable nuclei, which are on or near the /^-stability line 
in the nuclear chart. This is basically because only those 
nuclei have been accessible experimentally. In such sta- 
ble nuclei, the magic numbers suggested by Mayer and 
Jensen remain valid, and the shell structure can be under- 
stood well in terms of the harmonic oscillator potential 
with a spin-orbit splitting. 

Recently, studies on exotic nuclei far from the (3- 
stability line have started owing to the development of 
radioactive nuclear beams |2J. The magic numbers in 
such exotic nuclei can be a quite intriguing issue. In this 
Letter, we show that new magic numbers appear and 
some others disappear in moving from stable to exotic 
nuclei in a rather novel manner due to a particular part 
of the nucleon-nucleon interaction. 

In order to understand underlying single-particle prop- 
erties of a nucleus, we can make use of effective (spher- 
ical) single-particle energies (ESPE's), which represent 
mean effects from the other nucleons on a nucleon in 
a specified single-particle orbit. The two-body matrix 
element of the interaction depends on the angular mo- 
mentum J, coupled by two interacting nucleons in or- 
bits ji and ]2- Since we are investigating a mean effect, 
this J-dependence is averaged out with a weight factor 
(2 J +1), and only diagonal matrix elements are taken. 
Keeping the isospin dependence, T=0 or 1, the so-called 
monopole Hamiltonian is thus obtained with a matrix 
element 



(1) 



v t = Ej(2J + 1) < 3ij 2 \V\jiJ2 >JT 



Ej(2J+l) 



where < jiJ2\V\j' 1 j' 2 >jt stands for the matrix clement 



of a two-body interaction, V. 

The ESPE is evaluated from this monopole Hamilto- 
nian as a measure of mean effects from the other nucleons. 
The normal filling configuration is used. Note that, be- 
cause the J dependence is taken away, only the number 
of nucleons in each orbit matters. As a natural assump- 
tion, the possible lowest isospin coupling is assumed for 
protons and neutrons in the same orbit. The ESPE of an 
occupied orbit is defined to be the separation energy of 
this orbit with the opposite sign. Note that the separa- 
tion energy implies the minimum energy needed to take 
a nucleon out of this orbit. The ESPE of an unoccupied 
orbit is defined to be the binding energy gain by putting 
a proton or neutron into this orbit with the opposite sign. 

Figure shows neutron ESPE's for 30 Si and 24 0, both 
of which have A^=16. The Hamiltonian and the single- 
particle model space are the same as those used in [0, 
where the structure of exotic nuclei with N~20 has been 
successfully described within a single framework. The 
valence orbits are then 0^5/2,3/2, lsi/2, O/7/2 and Lj>3/2- 

The nucleus 30 Si has six valence protons in the sd shell 
and is a stable nucleus, while 24 has no valence proton 
and is a neutron-rich exotic nucleus. For 30 Si, the neu- 
tron 0^3/2 and lsi/2 are rather close to each other (see 
Fig. (a)). For 24 0, as shown in Fig. (b), the 0d 3/2 
is lying much higher and is quite close to the pf shell, 
giving rise to a large gap (~6 MeV) between 0d 3 / 2 and 
lsi/2 ||. On the other hand, for the stable nucleus 30 Si, a 
considerable gap (~4 MeV) is created between the 0(i 3 /2 
and the pf shell (See Fig. [j] (a)). Thus, the N=2Q magic 
structure is evident for 30 Si, whereas the 7V=16 magic 
number arises in 24 0. In 24 0, the 0d 3 / 2 is lying higher 
reflecting the large spin-orbit splitting which is basically 
the same as that for 17 0. Although this high-lying 0d 3 / 2 
orbit is not so relevant to the ground state of lighter O 
isotopes, it should affect binding energies of nuclei around 
24 0. Such an anomaly was pointed out by Ozawa et al. 
B in observed binding energy systematics. 

The dramatic change of ESPE's from 24 to 30 Si is pri- 



marily due to the strongly attractive interaction between 
a proton in 0d 5 / 2 and a neutron in 0d^/ 2 - A schematic 
picture on this point is shown in Fig. [j] (c) for the gen- 
eral cases comprised of a pair of orbits j> = 1+ 1/2 and 
j < = I — 1/2 with I being the orbital angular momentum. 
Note that j> and j < are nothing but spin-orbit coupling 
partners. The present case corresponds to 1—2. As one 
moves from 24 to 30 Si, six valence protons are put into 
the 0d 5 / 2 orbit. Consequently, due to the strong attrac- 
tion shown in Fig. YU (c), a neutron in 0<i 3 /2 is more bound 
in 30 Si, and its neutron 0d 3 / 2 ESPE becomes so low as 
compared to that in 24 where such attraction is absent. 
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FIG. 1. Neutron ESPE's for (a) 30 Si and (b) 24 0, relative 
to lsx/2- The dotted line connecting (a) and (b) is drawn to 
indicate the change of the 0d 3 / 2 level, (c) The major inter- 
action producing the basic change between (a) and (b). (d) 
The process relevant to the interaction in (c). 

The process illustrated in Fig. p] (d) produces the at- 
tractive interaction in Fig. [0 (c). The nucleon-nucleon 
(TV TV) interaction in this process is written as 



V Ta 



<rfr<r{r). 



(2) 



Here, the symbol "•" denotes a scalar product, r and 
a stand for isospin and spin operators, respectively, r 
implies the distance between two interacting nucleons, 
and f Ta is a function of r. In the long range (or no r- 
dependence) limit of f T a(r), the interaction in eq.(0) can 
couple only a pair of orbits with the same orbital angular 
momentum I, which are nothing but j > or j < . 

The a operator couples j> to j< (and vice versa) more 
strongly than j > to j > or j < to j < . Therefore, the spin 
flip process is more favored in the vertexes in Fig. [D (d), 
while spin non-flip contributions certainly exist. The 
same mathematical mechanism works for isospin: the r 
operator favors charge exchange processes. Combining 
these two properties, V T(T produces large matrix elements 
for the spin-flip isospin-flip processes: proton in j > — > 
neutron in j < and vice versa. This gives rise to the inter- 
action in Fig. pi (c) with a strongly attractive monopole 



term for the appropriate sign of V Ta . This feature is a 
general one and is maintained with /Vo-M in eq.(0) with 
reasonable r dependences. 

In stable nuclei with N~Z with ample occupancy of 
the j > orbit in the valence shell, the proton (neutron) j < 
orbit is lowered by neutrons (protons) in the j > orbit. In 
exotic nuclei, this lowering can be absent, and then the 
j < orbit is located rather high, not far from the upper 
shell. In this sense, the proton-neutron j>-j< interaction 
enlarges a gap between major shells for stable nuclei with 
proper occupancy of relevant orbits. This interaction has 
been known, for instance j9JJlO|], to play important roles 
also in other issues, e.g., the onset of the deformation. 

The origin of the strong V T<7 is quite clear. The One- 
Boson-Exchange-Potentials (OBEP) for 7r and p mesons 
have this type of terms as major contributions. While the 
OBEP is one of major parts of the effective NN interac- 
tion, the effective NN interaction in nuclei can be pro- 
vided by the G-matrix calculation with core polarization 
corrections and other various effects. Such effective NN 
interaction will be called simply G-matrix interaction for 
brevity. The G-matrix interaction should maintain the 
basic features of meson exchange processes, and, in fact, 
existing G-matrix interactions generally have quite large 
matrix elements for the cases shown in Fig. Pfl (c), includ- 
ing strongly attractive monopole terms 111]]. 

We would like to point out that the l/N c expansion of 
QCD by Kaplan and Manohar indicates that V T(J is one 
of three leading terms of the long-range part of the NN 
interaction [O. Since the next order of this expansion 
is smaller by a factor (l/iV c ) 2 , the leading terms should 
have rather distinct significance. One of the other two 



leading terms in the 1/N C expansion 12 is a central force, 



Vo = /o (r), 



(3) 



where /q is a function of r. The last leading term is a 
tensor force. We shall come back to these forces later. 

Figure || shows the effective 0d 3 / 2 -lsi/2 g a P, «-e., the 
difference between ESPE's of these orbits, in iV=16 iso- 
tones with Z=8~20 for three interactions: "Kuo" means 
a G-matrix interaction for the sd shell calculated by Kuo 
, and USD was obtained by adding empirical modifica- 
tions to "Kuo" [|J. The present shell- model interaction 
is denoted SDPF hereafter, and its sd-shell part is noth- 
ing but USD with small changes Q. Steep decrease of 
this gap is found in all cases, as Z departs from 8 to 
14. In other words, a magic structure can be formed 
around Z—8, but it should disappear quickly as Z de- 
viates from 8 because the gap decreases very fast. The 
slope of this sharp drop is determined by V Qd ~ ' od in 
eq. (Ty), where the dominant contribution is from T=Q. 
Note that Vl r ,= M 



is most attractive among all V/ ?2 's 



in the sd shell for each of the three interactions. In SDPF, 
for instance, its magnitude exceeds by ~1 MeV that of 
the second most attractive one. 



We shall now discuss briefly the present gap in other 
approaches. The gap can be calculated from the Woods- 
Saxon potential which is a standard modeling of single- 
particle structure for stable nuclei. The resultant gap 
is rather flat, and is about half of the SDPF value for 
Z=8. With Skyrme Hartree-Fock (HF) interactions, the 
gap changes more smoothly and gradually, too. 




FIG. 2. Effective lsi/2-0d 3 / 2 gap in iV=16 isotones as a 
function of Z. Shell model Hamiltonians, SDPF, USD and 
"Kuo" are used. See the text. 



The occupation number of the neutron lsx/% is cal- 
culated by Monte Carlo Shell Model [Q with full con- 
figuration mixing for the nuclei shown in Fig. 0. It is 
nearly two for 24 as expected for a magic nucleus, but 
decreases sharply as Z increases. It remains smaller (< 
1.5) in the middle region around Z=14, and finally goes 
up again for Z~20. This means that the N—1Q magic 
structure is broken in the middle region of the proton 
sd shell, where deformation effects also contribute to the 
breaking. The iV=16 magic number is thus quite valid 
at both ends. It is of interest that the gap becomes large 
again for larger Z, due to other monopole components. 

We now turn to exotic nuclei with 7V-20. The ESPE 
has been evaluated for them in Q. The small effective 
gap between Qd 3 / 2 and the pf shell for neutrons is ob- 
tained, and is found to play essential roles for various 
anomalous features. This small gap is nothing but what 
we have seen for 24 in Fig. [I] (b). Thus, the disappear- 
ance of N=2Q magic structure in Z=9~14 exotic nuclei 
and the appearance of the new magic structure in 24 
have the same origin. 

A very similar mechanism works for p-shell nuclei. We 
consider the structure of a stable nucleus 13 C, and exotic 
nuclei n Be and 9 He, all of which have N—7. Shell model 
calculations are performed with the PSDMK2 Hamilto- 
nian jl4] in the p+sd shell model space on top of the 4 He 
core. Figure indicates that the experimental levels are 
reproduced well for 13 C Jig ], whereas notable deviations 
are found in n Be and 9 He Kq,Eq|. Here, for experimen- 
tal levels in the continuum, their resonance(-like) energies 
are compared to the shell-model results. 

The p-shell part of PSDMK2 is the Cohen-Kurath 
(CK) Hamiltonian pTy. The single-particle energies of 



0p 3 / 2 and Opi/2 are 1.38 and 1.68 MeV, respectively, in 
PSDMK2. These energies correspond to the observed 
spectra of 5 He: 1.15 for 0p 3 / 2 and ~5 MeV for 0pi/ 2 . 
We use these observed values as single particle energies, 
while Opi/2 is quite low in PSDMK2. As compared to 
the G- matrix interaction uM, the interaction in Fig. p] 
(c) for the p-shell is too weak in CK, and is enlarged 
by shifting all < 0p 3/2 0pi/ 2 \V\0p 3/2 0p 1/2 >j,t=o's by 
-2 MeV independently of J. Thus, the Hamiltonian is 
modified only for three parameters. Since we are focus- 
ing on basic trends and would like to eliminate spurious 
center-of-mass components, all shell model calculations 
for N—7 isotones are made in the so-called Ohto or lhu> 
space. Figure || (a) indicates that the levels of 13 C cal- 
culated from the present Hamiltonian are similar to the 
ones obtained from PSDMK2, because the higher-lying 
Opi/2 is pulled down by V as discussed above. 
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FIG. 3. Energy levels for (a) 13 C, (a) n Be and (c) 9 He, 
relative to the state corresponding to the experimental ground 
state. The results by the PSDMK2 and present Hamiltonians 
are compared to experiments. Hatched area indicate contin- 
uum states with resonance(-like) structures. 

The ground state of n Be and 9 He is known for the in- 
version between 1/2+ and l/2~ [Q. The present calcu- 
lation reproduces this inversion for both nuclei, whereas 
the PSDMK2 fails in cither case. With the present 
Hamiltonian, neutron effective \s 1 1 2 - Qp\/ 2 gap decreases 
from 8.6 MeV for 13 C down to 0.8 MeV for 9 He. Dynam- 
ical correlations finally invert 1/2+ and l/2~ eigenstates 
p9[ . Thus, we achieve a reasonable description of stable 
and exotic nuclei with N=7. Although the result ap- 
pears promising, we have investigated only possible im- 
provements, and further studies are needed for an overall 
description of p-sd nuclei. For instance, in comparison to 



experiments, the present interaction gives a 3/2 state 
in 15 ^2 MeV too high fLq] , and 9 He is too unbound 
relative to 8 He by ~4.5 MeV @. The latter may be 
related to coupling to continuum and/or to nuclear-size 
dependence of two-body interaction. 

The neutron Opi/2 orbit becomes higher as the nucleus 
loses protons in its spin- flip partner 0p 3 / 2 . In nuclei such 
as He, Li and Be, the N=8 magic structure then disap- 
pears. In some cases, N—6 becomes magic: for instance, 
bound 8 He and unbound 9 He are obtained, similarly to 
bound 24 and unbound 25 0. 

As the neutron Opi/2 is shifted quite high in He, Li 
and Be, this orbit becomes close to the lsi/2 which may 
come down due to halo structure. This is a situation 
like the Efimov state POJ] . The pairing between these 
two orbits, including contribution from the tensor force, 
may stabilize the two-neutron halo and provide us with 
bound n Li, etc. Thus, the present mechanism plays a 
crucial role for the structure of dripline nuclei. Without 
this mechanism, it is usually difficult to make Opi/2 and 
Is 1/2 so close to each other, except for an empirical shift 
of the single-particle energy. 

Moving back to heavier nuclei, from the strong interac- 
tion in Fig. p] (c), we can predict other magic numbers, for 
instance, iV=34 associated with the O/7/2 - Of 5/2 interac- 
tion. In heavier nuclei, 0(77/2, Oft.9/2, etc. are shifted up- 
ward in neutron-rich exotic nuclei, disturbing the magic 
numbers N—82, 126, etc. It is of interest how the re- 
process of nucleosynthesis is affected by it. 

In conclusion, we showed how magic numbers are 
changed in nuclei far from the /3-stability line: N—6, 16, 
34, etc. can become magic numbers in neutron-rich ex- 
otic nuclei, while usual magic numbers, N=8, 20, 40, 
etc., may disappear. Since such changes occur as re- 
sults of the nuclear force, there is isospin symmetry that 
similar changes occur for the same Z values in mirror 
nuclei. The mechanism of this change can be explained 
by the strong V Ta interaction which has robust origins 
in OBEP, G-matrix and QCD. In fact, simple structure 
such as magic numbers should have a simple and sound 
basis. Since it is unlikely that a mean central potential 
can simulate most effects of V Ta , we should treat V Ta 
rather explicitly. It is nice to build a bridge between 
very basic feature of exotic nuclei and the basic theory of 
hadrons, QCD. In existing Skyrme HF calculations ex- 
cept for those with Gogny force, effects of V Ta may not 
be well enough included, because the interaction is trun- 
cated to be of (5-function type. The Relativistic Mean 
Field calculations must include pion degrees of freedom 
to be consistent with V Ta . Thus, the importance of V Ta 
opens new directions for mean field theories of nuclei. In 
addition, fo{f) in eq.(0) and f T(T {r) in eq.(H) should have 
different ranges. Since the latter has smaller contribu- 
tions in exotic nuclei, this difference should produce very 
interesting effects on nuclear size. Loose-binding or con- 
tinuum effects are important in some exotic nuclei. By 



combining such effects with those discussed in this Letter 
one may draw a more complete picture for the structure 
of exotic nuclei. Finally, we would like to mention once 
more that the V Ta interaction should produce large, sim- 
ple and robust effects on various properties, and may 
change the landscape of nuclei far from the /^-stability 
line in the nuclear chart. 
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